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ABSTRACT 

m 

\ We have searched for molecular hydrogen in damped Lyman-a (DLA) and sub- 

■ DLA systems at high redshift (zabs > 1-8) using UVES at the VLT down to a detection 

' limit of typically A^(H2) = 2 x lO^"* cm~^. Out of the 33 systems in our sample, 8 have 

firm and 2 have tentative detections of associated H2 absorption lines. Considering that 
^— ^ ' 3 detections were already known from past searches, molecular hydrogen is detected in 

^r~i \ 13 to 20 percent of the newly-surveyed systems. We report new detections of molecular 

(3 ! hydrogen at ^abs = 2.087 and 2.595 toward, respectively, Q 1444+014 and Q 0405-443, 

and also reanalyse the system at Zabs = 3.025 toward Q 0347— 383. 

In all of the systems, we measure metallicities relative to Solar, [X/H] (with either 
X=Zn, or S, or Si), and depletion factors of Fe, [X/Fe], supposedly onto dust grains, 
O ' and compare the characteristics of our sample with those of the global population 

of DLA systems (60 systems in total). We find that there is a correlation between 
c/3 ■ metallicity and depletion factor in both our sample and also the global population of 

DLA systems. Although H2 molecules are detected in systems with [Zn/Fe] as small 
as 0.3, the DLA and sub-DLA systems where H2 is detected are usually amongst 
those having the highest metallicities and the largest depletion factors. In particular, 
H2 is detected in the five systems having the largest depletion factors. Moreover, 
i the individual components where II2 is detected have depletion factors systematically 

larger than other components in the profiles. In two different systems, one of the H2- 
detected components even has [Zn/Fe] ^ 1.4. These are the largest depletion factors 
ever seen in DLA systems. All this clearly demonstrates the presence of dust in a large 
fraction of the DLA systems. 

The mean H2 molecular fraction, / = 27V(H2)/[27V(H2) + A^(Hi)], is generally 
small in DLA systems (typically log/ < —1) and similar to what is observed in the 
Magellanic Clouds. There is no correlation between the observed amount of II2 and the 
Hi column density. In fact, two systems where II2 is detected have log (Hi) < 20.3 
and, therefore, are sub-DLA systems. From 58 to 75 percent of the DLA systems have 
log/ < —6. This can be explained if the formation rate of H2 onto dust grains is 
reduced in those systems, probably because the gas is warm {T > 1000 K) and/or the 
ionizing flux is enhanced relative to what is observed in our Galaxy. 

Key words: Cosmology: observations - Galaxies: haloes - Galaxies: ISM - Quasars: 
absorption lines - Quasars: individual: Q 0347-383, Q 0405-443, Q 1444-^014 



o 

> 



* Based on observations carried out at the European Southern 
Observatory (ESO), under visitor mode progs. ID 65.O-0063, 
66.A-0624, 67.A-0078, 68.A-0106 and 68.A-0600, with the UVES 



echelle spectrograph installed at the ESO Very Large Telescope 
(VLT), unit Kuoyon, on mount Paranal in Chile. 
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1 INTRODUCTION 

High-redshift damped Lyman-a (DLA) systems observed in 
absorption against QSO spectra are characterized by their 
extremely strong Hi A1215 lines corresponding to large neu- 
tral hydrogen column densities, A''(Hl) ^ 2x 10^" cm"'^. Hy- 
drodynamical simulations suggest that DLA systems are lo- 
cated inside regions of over-densities of the order of 1000 and 
higher (see e.g. Haehnelt, Steinmetz & Ranch 2000; Gardner 
et al. 2001) and that DLA systems at Zabs Si 2 occur very 
close (within 10-15 kpc) to the center of L* galaxies. This 
is purely speculative however and the exact nature of DLA 
systems is still to be clarified. Though observational studies 
of DLA systems have been pursued over two decades now, 
important questions are still unanswered, such as: (i) the 
presence of in-situ star-formation activity in DLA systems, 
(ii) the connection between observed abundance ratios and 
dust content, (iii) how severe is the bias due to dust ob- 
scuration in current DLA samples. Assessing the molecular 
content of DLA systems can provide direct handle on at 
least some of these issues. 

Formation of H2 is expected on the surface of dust 
grains if the gas is cool, dense and mostly neutral, and from 
the formation of negative hydrogen if the gas is warm and 
dust-free (see e.g. Jenkins & Peimbert 1997; Cazaux & Tie- 
lens 2002). As the former process is most likely dominant 
in the neutral gas associated with DLA systems, it is possi- 
ble to obtain an indirect indication of the amount of dust in 
DLA systems without depending on extinction and/or heavy 
element dust depletion effects. Moreover, from determining 
the populations of different rotational levels of H2, one can 
constrain kinetic and rotational excitation temperatures as 
well as particle densities. The effective photo-dissociation of 
H2 takes place in the energy range 11.1 — 13.6 eV through 
Lyman- and Werner-band absorption lines and, therefore, 
the intensity of the local UV radiation field can be derived 
from the observed molecular fraction. A direct determina- 
tion of the local UV radiation field could have important 
implications in bridging the link between DLA systems and 
star-formation activity in high-redshift galaxies. 

At Hi column densities as large as those measured in 
DLA systems, H2 molecules are conspicuous in our Galaxy: 
gas clouds with log A''(Hi) > 21 usually have log A''(H2) > 19 
(see e.g. Savage et al. 1977; Jenkins & Shaya 1979). Given 
this fact, it is somewhat surprising that earlier searches for 
molecular hydrogen in DLA systems, though not systematic, 
have found either small values or upper limits on the molec- 
ular fraction of the gas, / = 2iV(H2)/[2iV(H2) + iV(Hl)] 
(Black, Chaffee & Foltz 1987; Levshakov et al. 1992). For 
a long time, only the DLA system at z^ha = 2.811 to- 
ward Q 0528—250 was known to contain molecular hydro- 
gen (Levshakov & Varshalovich 1985; Foltz, Chaffee & Black 
1988). More recently, Ge & Bechtold (1999) searched for H2 
molecules in a sample of eight DLA systems using the MMT 
moderate-resolution spectrograph (FWHM = 1 A). Apart 
from the detection of molecular hydrogen at Zabs = 1.973 
and 2.338 toward, respectively, Q 0013-004 and Q 1232+082 
(Ge & Bechtold 1997; Ge, Bechtold & Kulkarni 2001), they 
measured in the other systems upper limits on / in the range 
10"^ - 10"*. 

The molecular hydrogen content of the above three 
systems, where H2 was detected at intermediate spectral 



resolution, has been reexamined systematically using high 
spectral resolution data (FWHM = 0.1 A), leading to: 
iV(H2) ~ 6x10^® cm"^ and / ~ 5x10"^ toward Q 0528-250 
(Srianand & Petitjean 1998), 5 x 10^^ < iV(H2) < 10^" 
cm"^ and 2 x 10"^ < f < 0.2 toward Q 0013-004 (Pe- 
titjean, Srianand & Ledoux 2002) and A''(H2) > 2 x 10^'' 
cm"2 and / > 4 X 10"" toward Q 1232-1-082 (Srianand, Pe- 
titjean & Ledoux 2000). Recent new detections are from 
Levshakov et al. (2002) at Zabs = 3.025 toward Q 0347-383 
and Ledoux, Srianand & Petitjean (2002b) at Zabs = 1.962 
toward Q 0551—366. Two tentative detections have also been 
reported at Zabs = 2.374 toward Q084H-129 (Petitjean, 
Srianand & Ledoux 2000) and at Zabs = 3.390 toward 
Q 0000—263 (Levshakov et al. 2000) based, however, on the 
identification of two weak features located in the Lyman-a 
forest. 

This paper describes the high spectral resolution sur- 
vey of high-redshift DLA systems, focused on the search for 
molecular hydrogen, that we have carried out at the ESO 
VLT during the first two years of operations of the UVES 
echelle spectrograph. We present here our final results from 
the observations of a total sample of 24 DLA (logA'^(Hl) ^ 
20.3) and 9 sub-DLA (19.4 < logAf(Hl) < 20.3) systems, 
including two new H2 detections. Preliminary results from 
the observations of a sample of 11 systems were presented 
in Petitjean et al. (2000). Sect. |5|in the present paper de- 
scribes the observations, the data and the absorber sample, 
Sect.|3gives details on three DLA/sub-DLA systems where 
H2 is detected, and in both Sects. |l]and|3general results are 
discussed from the study of the whole sample. We conclude 
in Sect.|S| 



2 OBSERVATIONS AND SAMPLE 

The Ultraviolet and Visible Echelle Spectrograph (UVES; 
Dekker et al. 2000), installed at the ESO VLT 8.2-m tele- 
scope, unit Kueyen, on Mount Paranal in Chile, was used 
during five visitor-mode runs, from April 2000 to January 
2002. During each run, we surveyed DLA and sub-DLA sys- 
tems with the aim of searching for H2 lines with redshifts in 
the range 1.8 — 3.4. In a first step, we observed the {V < 19) 
background QSOs twice for 1 to 1.5 hr over wavelength 
ranges covering the location of the main metal lines of the 
absorption systems and possibly associated H2 features in, 
respectively, the Red and the Blue arms of UVES (which 
were used simultaneously). For most of the targets, this re- 
sulted to a signal-to-noise ratio in excess of 10 in the Blue. 
The quick-look data reduction package available at the tele- 
scope was extensively used to decide in real time whether H2 
absorption lines are present or not. This allowed us to opti- 
mize our observing strategy: in case of detection, additional 
exposures were gathered to achieve the high signal-to-noise 
ratio needed for an accurate determination of the H2 column 
densities in different rotational levels. In addition to our ob- 
servations, we also included the few UVES archival spectra 
of QSOs with DLA systems for which the H2 wavelength 
ranges were covered. 

The total QSO sample presently under consideration is 
made of 25 lines of sight. Observations toward Q 0013—004, 
Q 0551-366 and Q 1232-1-082, along the lines of sight of 
which H2 is detected, are described respectively in Petit- 
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Table 1. Atomic data 



Transition 


Avac 

(A) 


/ 


Hei. 


Transition 


Avac 

(A) 


/ 


Ker. 


H I A1025 


1025.7223 


0.07912 


a 


Sin A1808 


1808.0129 


0.00208 


f 


H1AI215 


1215.6701 


0.4164 


a 


P II A963 


963.801 


1.458 


a 


Ci A1328 


1328.8333 


0.0630 


b 


P11AII52 


1152.8180 


0.236 


a 


Ci*A1329.08 


1329.0853 


0.0213 


b 


Sii A1250 


1250.584 


0.00545 


a 


Ci*A1329.10 


1329.1004 


0.0260 


b 


Sii A1253 


1253.811 


0.0109 


a 


Ci*A1329.12 


1329.1233 


0.0160 


b 


Sii A1259 


1259.519 


0.0162 


a 


CiA1560 


1560.3092 


0.0719 


b 


Ar I A1048 


1048.2199 


0.257 


g 


Ci*A1560.6 


1560.6822 


0.0539 


b 


Tin A1910.6 


1910.609 


0.104 


h 


Ci*A1560.7 


1560.7090 


0.0180 


b 


Tin A1910.9 


1910.938 


0.098 


h 


Ci A1656 


1656.9283 


0.139 


b 


Cr II A2056 


2056.2569 


0.105 




Ci*A1656 


1656.2672 


0.0589 


b 


Cr II A2062 


2062.2361 


0.0780 




Ci*A1657.3 


1657.3792 


0.0356 


b 


Cr II A2066 


2066.1640 


0.0515 




Ci*A1657.9 


1657.9068 


0.0473 


b 


Fe II A1096 


1096.8769 


0.0324 




Ci**A1657 


1657.0082 


0.104 


b 


Fell A1121 


1121.9749 


0.0202 




Cii*A1037 


1037.0182 


0.123 


a 


Fell A1125 


1125.4478 


0.0163 




Cii*A1335.6 


1335.6627 


0.0128 


a 


Fell A1143 


1143.2260 


0.0177 




Cii*A1335.7 


1335.7077 


0.115 


a 


Fell Al 144 


1144.9379 


0.106 




N1AII34.I 


1134.1653 


0.0134 


a 


Fe II A1608 


1608.4509 


0.0585 




N1AII34.4 


1134.4149 


0.0268 


a 


Fell A1611 


1611.2005 


0.00136 




N1AII34.9 


1134.9803 


0.0402 


a 


Fe II A2249 


2249.8768 


0.00182 




N1AI2OO.2 


1200.2233 


0.0885 


a 


Fe II A2260 


2260.7805 


0.00244 




1 A950 


950.8846 


0.00157 


a 


Fe II A2374 


2374.4612 


0.0313 




1 A974 


974.07(5) 1 


0.00002 


a 


Ni IIA1317 


1317.2170 


0.0774 


k 


MgiiA1239 


1239.9253 


0.000554 


c 


Ni II A1370 


1370.1320 


0.0765 


k 


MgiiA1240 


1240.3947 


0.000277 


c 


NiiiA1741 


1741.5531 


0.0427 


1 


Alii A1670 


1670.7886 


1.833 


a 


NiiiA1751 


1751.9157 


0.0277 


1 


Si 11AIO2O 


1020.6989 


0.0283 


a 


Znll A2026 


2026.1371 


0.489 


i 


Si IIA1304 


1304.3702 


0.0894 


d 


Zn II A2062 


2062.6604 


0.256 


i 


Si IIA1526 


1526.7070 


0.127 


e 











1 See Subsect. |0] 

2 References for oscillator strengths: (a) Morton (1991); (b) Wiese, Fuhr & Deters (1996); (c) Welty et al. (1999); (d) Spitzer & 
Fitzpatrick (1993); (e) Scliectman, Povolny & Curtis (1998); (f) Bergeson & Lawler (1993b); (g) Federman et al. (1992); (h) Wiese, 
Fedchak & Lawler (2001); (i) Bergeson & Lawler (1993a); (j) Howk et al. (2000); (k) Fedchak & Lawler (1999); (1) Fedchak, Wiese & 
Lawler (2000). 



jean et al. (2002), Ledoux et aL (2002b) and Srianand 
et al. (2000). In addition, new data for Q 0528-250 and 
Q 1232+082 were acquired recently and the corresponding 
H2-detected systems should be analysed in detail in forth- 
coming papers. However, we give in Sect.|l]the measurement 
from UVES data of the molecular hydrogen and metal con- 
tents of the DLA system at ^abs = 2.811 toward Q 0528—250. 

With respect to the remaining and/or new detections 
of H2 in our sample, high-resolution, high signal-to-noise 
ratio spectra of Q 0347— 383 (zcm ~ 3.21, mv = 17.8), 
Q 0405-443 (zom = 3.02, tub = 17.6) and Q 1444+014 
(zom = 2.21, rriB, = 18.5) were obtained on, respectively, 
January 7-9, 2002, October 20-23, 2001, and June 16, 2001. 
During the observations of Q 0405—443, central wavelengths 
3900 and 5700 A were used in, respectively, the Blue and 
the Red arms of UVES with Dichroic #1, and central wave- 
lengths were adjusted to 4370 A in the Blue and 7400 A 
in the Red with Dichroic #2. Full wavelength coverage was 
achieved this way from 3283 to 7320 A and from 7468 to 9302 
A accounting for the gap between the two Red-arm CCDs. 
During the observations of Q 1444+014 (resp. Q 0347-383), 
central wavelengths were adjusted to 3800 A (resp. 4300 A) 
in the Blue and 5640 A (resp. 8000 A) in the Red, and we 



used Dichroic #1 (resp. #2). Complementary UVES data 
for Q 1444+014, obtained on May 29-30, 2000, with sim- 
ilar settings, were retrieved from the UVES archive. For 
each of the 25 quasars in our sample, the CCD pixels were 
binned 2 x 2 in both arms and the slit width was fixed to 1" 
(except for Q 0347— 383: 0.8") yielding a spectral resolution 
R = 42500 (resp. R = 52000). The total integration times 
were 9 x 4500 s for Q 0347-383, 4 x 4350 s with Dichroic 
#1 (and 2 x 4350 s with Dichroic #2) for Q 0405-443, and 
8 x 3600 s for Q 1444+014. 

The spectra for each QSO were reduced in the dedicated 
context of MIDAS, the ESQ data reduction system, using 
the UVES pipeline (Ballester et al. 2000) in an interactive 
mode. The main characteristics of the pipeline are to per- 
form a precise inter-order background subtraction for science 
frames and master flat-fields, and an optimal extraction of 
the object signal rejecting cosmic ray impacts and subtract- 
ing the sky spectrum simultaneously. The pipeline products 
were checked step by step. The wavelength scale of the spec- 
tra reduced by the pipeline was then converted to vacuum- 
heliocentric values and individual 1-D exposures were scaled, 
weighted and combined altogether using the scombine task 
of the NOAO onedspec package of IRAF. During this pro- 
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Figure 1. Velocity profiles of lines from low ions and neutral species in the DLA system at z^bs = 3.025 toward Q 0347—383. The 
Sll A1259 profile is reproduced from Prochaska & Wolfe (1999). Our best-fitting model is superimposed on the spectra with vertical lines 
marking the location of individual components. The dotted parts in some of the synthetic profiles correspond to other transitions than 
the ones indicated (ZnnA2062 and CrllA2062). The single component where H2 is detected is labelled with number 1. The location of 
telluric absorption features is marked by the symbol ©. Note that the blending with telluric lines has been carefully taken into account 
(dark curves) for the fitting of the Sin A1808 and Znll A2026 profiles (light curves). 
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Figure 2. Voigt-profile fitting to the transition lines from the J = 0, 1, 2, 3 and 4 rotational levels of the vibrational ground-state Lyman 
and Werner bands of H2 at z^bs = 3.02489 toward Q 0347— 383. Only a few of the observed bands are displayed. H2 is detected on this 
line of sight in a single gas cloud. 



cess, the spectra of each QSO were rebinned to a common 
wavelength step (0.05 A pix~^ or less). In particular, the 
spectra of Q 0347-383, Q 0405-443 and Q 1444+014, where 
H2 is detected, were rebinned to, respectively, 0.03, 0.0471 
and 0.0339 A pix ^ , yielding after combination of individual 
exposures signal-to-noise ratios in, respectively, the ranges 
20-50, 10-50 and 10-45 in the Blue and 90-15, 120-20 and 
90-30 in the Red. 

Standard Voigt-profile fitting techniques were used to 
determine column densities, in all the absorption line sys- 
tems for which the H2 wavelength range was observed, using 
the oscillator strengths compiled in Tabled for metal ions 
and the oscillator strengths given by Morton & Dinerstein 
(1976) for H2. In the following, we measure metallicities rela- 
tive to Solar, [X/H] = log[iV(X)/iV(H)] - log[iV(X)/iV(H)]0, 
with either X=Zn, or S, or Si, and depletion factors of iron, 
[X/Fe] = log[Af(X)/iV(Fe)] - log[iV(X)/Af(Fe)]0, adopting 
the Solar abundances from Savage & Sembach (1996). 

There has been no a priori selection of the systems. 
In our sample, the 24 DLA and 9 sub-DLA absorbers are 
generally drawn at random from the global population of 
known systems. However, we have reobserved with UVES 
three DLA systems where II2 was previously detected (to- 
ward Q 0013-004, Q 0528-250 and Q1232-f082). Finally, 



the number of lines of sight along which several systems are 
seen is larger in our sample than in the global population 
of DLA systems. One would expect this to introduce a bias 
in favor of low-metallicity absorbers. However, the charac- 
teristics of our sample and those of the global population of 
DLA systems are compared in Section |1| and they are shown 
to be similar. 



3 DESCRIPTION OF INDIVIDUAL SYSTEMS 
WITH H2 DETECTION 

In this Section, we reconsider the case of the DLA sys- 
tem at 2;abs = 3.025 toward Q 0347—383 and describe 
two new detections of molecular hydrogen at, respectively, 
Zabs = 2.595 and 2.087 toward, respectively, Q 0405-443 
and Q 1444-^014. 

3.1 Q 0347-383, z^bs = 3.025 

The detection of molecular hydrogen toward this quasar has 
been reported recently by Levshakov et al. (2002) who used 
UVES commissioning data. As part of our systematic survey, 
we have obtained significantly better UVES data for this line 
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Table 2. Voigt-profilc fitting results for different rotational levels 
of tlie vibrational ground-state Lyman and Werner bands of H2 
toward Q 0347-383 



Table 3. Ionic column densities in individual components of the 
DLA system at z^bs = 3.025 toward Q 0347-383 



■^abs 


Level 


log N ± (Tiog jv 


b±o-j 






(H2) 


(km s"-*-) 


3.02489 


J = 


13.25 ±0.08 


1.3 ±0.2 




J = 1 


14.26 ±0.06 






J = 2 


13.65 ±0.04 






J = 3 


13.90 ±0.04 






J = 4 


13.12 ±0.12 






J = 5 


^ 12.75 





Possible blends. 

of sight and we reanalyse here the DLA system at ^abs = 
3.025 in the same footing as the rest of our sample. 

We display in Fig. Q the velocity profiles of the most 
important lines from low ions and neutral species that are 
present in our data. Two well-separated components are ob- 
served for a number of ions, namely Nl, Ol, Sill, Pll, Sll 
and An. One of these two components, at Zabs ~ 3.02485, 
has detected H2 lines. 

Fitted altogether (see below), the detected H2 lines are 
found to be consistently located at a redshift of 3.02489 
and show a systematic offset of ~ 3 km s~^ relative to 
the deepest part of the low-ion profiles. We simultaneously 
fitted the lines from different rotational levels using the 
same broadening parameter and redshift, and measured ac- 
curate column densities for the J = to 4 levels from sub- 
samples of 4 to 7 unblended lines per J level. The results 
are given in Table Inland the best-fitting model is shown in 
Fig. |5| Note in particular that the measurement of the col- 
umn density in J = is based on the detection at more 
than 5a of four transition lines: Wl-OR(O), L7-0R(0), L3- 
0R(0) and L2-0R(0). A 3a absorption feature, possibly 
corresponding to C I A1656, is also detected at a redshift 
consistent with that of the lines. If this is true then 
logAf(Cl) = 11.73 ± 0.26 (see Table EJ. The second low- 
ion line component, at Zabs = 3.02463, has no detected 
H2 down to, respectively, logAf(H2) < 13.0 and < 13.4 
(5a limits) for, respectively, J = and 1. We measured in 
this DLA system a total neutral hydrogen column density 
log iV(H l) = 20.56 ± 0.05 (see Fig. and a total molec- 
ular hydrogen column density logA^(H2) — 14.55 ± 0.09 
(see Table This leads to an overall molecular fraction, 
log/ = —5.71 ± 0.10, in agreement with that derived by 
Levshakov et al. (2002). 

We measured ionic column densities in the two main 
components of the system and in two weaker ones only for 
species with unsaturated and/or mildly saturated lines. The 
results are given in Table|3 The broadening parameters were 
mainly constrained by Fen and, to some extent also, Nl 
lines covering a range of oscillator strengths. It can be seen 
in Fig. that in our data the Znll doublet lines are firmly 
detected in between sky lines in the strongest component of 
the system (zabs = 3.02485) and are also likely present in the 
component at Zabs ~ 3.02463. In addition, the Crll triplet 
is nicely detected in each of the four considered components 
(see Fig. [Hand Table gj. 

In Fig. we show evidence that the SillA1808 line 
is badly blended with a telluric absorption feature which 



Ion 


Transition 


log N ± (Tiog JV 


b±at 




lines used 




(km s"-*-) 




3.02434 






Ni 


1134.4,1134.9,1200.2 


13.14 ±0.08 


6.7 ±0.8 


Cr 11 


2056,2062,2066 


12.09 ±0.14 




Fell 


1096,1125,1144 


13.36 ±0.06 




Nin 


1741,1751 


12.37 ± 0.36 


„ 


^abs 


3.02463 






Nl 


1134.4,1134.9,1200.2 


14.15 ±0.03 


11.4 ±0.8 


Oi 


950,974 


16.12 ±0.12 ^ 




Sin 


1020,1808 


14.47 ±0.03 




Pii 


963 


12.34 ±0.19 " 




Sii 


1259 


14.50 ±0.03 




An 


1048 


13.68 ±0.03 




Ti 11 


1910.6,1910.9 


< 11.86 <= 




Cr 11 


2056,2062,2066 


12.69 ±0.05 




Fen 


1096,1125,1144 


14.20 ±0.02 




Nin 


1741,1751 


13.07 ± 0.09 




Znn 


2026,2062 


11.81 ±0.12 




^abs — 








Ci 


1656 ' 


11 7"? + 26 




Ci* 


1656 


< 11.75 " 




Ci** 


1657 


< 11.50 ° 




Cii* 


1037 


13.55 ± 0.23 


4.9 ± 0.3 


Nl 


1134.4,1134.9,1200.2 


14.47 ± 0.03 




Oi 


950,974 


16.18 ±0.18 ^ 




Si II 


1020,1808 


14.48 ± 0.04 




Pii 


963 


12.48 ± 0.09 




Sii 


1259 


14.36 ± 0.04 




Ar I 


1048 


13.44 ± 0.05 


„ 


Ti II 


1910.6,1910.9 


< 11.86 <= 




Crn 


2056,2062,2066 


12.29 ±0.11 




Fen 


1096,1125,1144 


13.81 ±0.05 




Nin 


1741,1751 


12.65 ±0.19 




Znn 


2026,2062 


12.02 ±0.04 




■^abs 


3.02501 






Ni 


1134.4,1134.9,1200.2 


13.44 ±0.04 


6.1 ± 0.5 


Sii 


1259 


13.51 ±0.12 




An 


1048 


12.71 ±0.06 




Cr II 


2056,2062,2066 


12.25 ±0.09 




Fell 


1096,1125,1144 


13.69 ±0.03 




Ni II 


1741,1751 


12.56 ±0.17 





Also 1608, 1611. 

An additional error of 20 percent was incorporated to reflect 
the uncertainty on the oscillator strength of the O I A974 line 
(see Morton 1991). 

The quoted error takes into account the main uncertainty 
which is related to the continuum placement because the single 
line used is located in a crowded Lyman-o forest. 

Line observed by Prochaska & Wolfe (1999). 
" 5cr upper limit. 

' A 3cr absorption feature is detected at Zabs = 3.02490 (~ 3 km 
s~^ redward of ^abs = 3.02485) while the H2 lines are 
consistently fitted with a single component at ^abs = 3.02489 
(see Table 13 . 
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Figure 3. Portions of the normalized UVES spectrum showing 
the damped Lya and Ly/3 lines of the DLA system at z^^i^s = 3.025 
toward Q 0347— 383. The best-fitted Voigt profiles superimposed 
on the data correspond to A''(Hl) = (3.6 ±0.4) X 10^° cm"^. 



we have carefully removed via modelling. The best-fitting 
of the SnA1259 line from the Keck-HIRES spectrum of 
Prochaska & Wolfe (1999; private communication) and of 
both SillAA1020,1808 lines in our VLT-UVES spectrum 
implies that the a- over iron-peak elemental ratio derived 
from the two main components of the system is close to 
Solar and may even be slightly under-Solar: [S, Si/Zn] — 
-0.11 ± 0.09, -0.35 ± 0.09. This is line with the findings 
of studies of large DLA samples (Prochaska & Wolfe 2002; 
Ledoux, Bergeron & Petitjean 2002a). In addition, we re- 
port the first detection of the O I A974 line in a DLA system 
(at the 4(T significance level; see Fig.0. The vacuum wave- 
length and the oscillator strength of this line are poorly- 
known quantities however (see Morton 1991). In order to fit 
the line, we used Avac = 974.075 A, and we included at the 
end an extra 20 percent uncertainty on the measured O I 
column density to reflect the uncertainty on the oscillator 
strength. We obtain [O/Zn] = 0.00 ± 0.20. Being Solar, this 
ratio is in agreement with the ones previously derived for S 
and Si. This is at variance with the claim by Levshakov et 
al. (2002) that in this DLA system a-elements are strongly 
overabundant compared to zinc ([a/Zn] — 0.6 ± 0.1). The 
latter authors have underestimated both the contamination 
of metal lines by telluric absorption features at Aobs > 7000 
A and also the saturation of O I lines in the Lyman-a forest. 

We compute metallicities [Zn/H] = -0.98 ± 0.09 and 
[S/H] = -1.09±0.06, and a mean depletion factor [Zn/Fe] = 
0.74 ± 0.09, using the integrated column densities from the 
two main components of the system at Zabs = 3.02463 and 
3.02485. Interestingly, the depletion, supposedly onto dust 
grains, and maybe also the metallicity, is larger in the com- 
ponent where H2 is detected ([Zn, S/Fe] = 1.07±0.06, 0.79± 
0.06 at ^abs = 3.02485) relative to the other components 
([Zn, S/Fe] = 0.47±0.12, 0.54±0.04 and [S/Fe] = 0.06±0.12 



at, respectively, Zabs = 3.02463 and 3.02501). Fig. ^indeed 
strikingly shows that the absorption line profiles of little or 
non-depleted elements, such as Si, P and S, are different 
from, e.g., the Cr and Fe profiles in their being about twice 
as narrow. 

3.2 Q 0405-443, Zabs = 2.595 

Lopez et al. (2001) reported the discovery of three DLA 
candidates in a low-resolution spectrum of Q 0405—443. We 
confirm here from UVES spectroscopy the damped nature 
of these absorption systems at redshifts Zaba = 2.550, 2.595 
and 2.621. Because the continuum of this quasar is especially 
well defined in our high-resolution spectra, we are able to 
derive accurate total neutral hydrogen column densities for 
each absorber. For this, we fitted the Lyman series from Lya 
to Ly6 and obtain logiV(Hl) = 21.00 ± 0.15, 20.90 ± 0.10 
and 20.25 ± 0.10 at, respectively, the above redshifts. Fig. |^ 
shows the best-fitting to the data which is in practice mainly 
constrained by the combination of Lya and Ly/9. A complete 
analysis of the three absorption systems will be published 
in Lopez et al. (in prep.) who have independently acquired 
UVES data of a similar quality for the purpose of studying 
elemental abundances. We concentrate here on only some of 
the metal lines observed in the DLA system at z^hs = 2.595 
(see Fig. 01 as in this particular absorber we found H2 . 

Molecular hydrogen at Zabs ~ 2.595 is detected in the 
J = 0, 1, 2 and 3 rotational levels. Interestingly, quite a 
large number of II2 lines from the L = to 14 Lyman-bands 
is observed at signal-to-noise ratios of the order of 10 to 30 
from Aobs = 3500 to 4000 A. This leads to a precise determi- 
nation of the molecular hydrogen column densities in J = 
and 1, and to some extent also, in J = 2 and 3 (see below 
however). Unsaturated and/or only moderately strong lines 
are observed for all J ^ 3 and damping wings are present in 
Werner-band lines from J — and 1. We simultaneously fit- 
ted all unblended lines from different rotational levels, using 
the same broadening parameter and redshift, with a single 
component at Zabs = 2.59471. The column density in each J 
level was derived from several trials of Voigt-profile fitting 
to take into account the range of possible b values. Errors 
in the column densities therefore correspond to a range of 
column densities and not to the rms error from fitting the 
Voigt profiles. The results are given in Table0|and the best- 
fitting model is shown in Fig. |S] The broadening parameter 
is found to be 6 = 1.5 km s~^. We note however that the 
relative strengths of the J ^ 2 lines require 6 > 1.5 km 
s~^ whereas the J ^ 1 lines are consistent with b < 1.5 km 
s~^. The broadening parameter could thus be smaller for 
lines from lower J and larger for lines from higher J, which 
may affects the determination of the column densities only 
in J ^ 2. 

The total molecular hydrogen column density measured 
in this system, logA''(H2) = 18.16lQgg, is one of the largest 
ever seen in DLA systems (see Sect. [IJ. However, the cor- 
responding molecular fraction, log / = — 2.44]IiQ'^2 1 rela- 
tively small due to the fairly large total neutral hydrogen 
column density. In spite of the high signal-to-noise ratio in 
our data at the location of the C l A1656 transition line, neu- 
tral carbon is not clearly detected. A > 5cr absorption fea- 
ture, whose origin is uncertain, is present at Zabs ~ 2.59485, 
corresponding to logA^(Cl) ^ 11.95. 
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Figure 4. Portions of the normalized UVES spectrum showing the damped Lya and Ly/3 lines of the DLA/sub-DLA systems at 
2abs = 2.550, 2.595 and 2.621 (see vertical lines) toward Q 0405— 443. The best-fitted Voigt profiles superimposed on the data correspond 
to NCR i) = (10 ± 4, 7.9 ± 2.1, 1.8 ± 0.5) x 10^0 cm-^ for each of the systems respectively. 
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Figure 5. Velocity profiles of some of the low-ion lines in the DLA system at ^abs = 2.595 toward Q 0405— 443. Our best-fitting model 
is superimposed on the spectra with vertical lines marking the location of individual components. The single component where H2 is 
detected is labelled with number 1. This component is clearly seen in the Cll*A1335.7 transition line and, to some extent also, in the 
Sll and Znll profiles. 
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Figure 6. Voigt-profile fitting to the transition lines from the J = 0, 1, 2 and 3 rotational levels of the vibrational ground-state Lyman 
and Werner bands of H2 at 2abs = 2.59471 toward Q 0405— 443. Only a few of the observed bands are displayed. H2 is detected on this 
line of sight in a single gas cloud. 



Table 4. Voigt-profile fitting results for different rotational levels 
of the vibrational ground-state Lyman and Werner bands of H2 
toward Q 0405-443 



■^abs 


Level log A'' di cTiog AT 






(H2) 


(km s~^) 



2.59471 



J = 

J = 1 

J = 2 

J = 3 

J = 4 

J = 5 



17.73 



-1-0.15 
0.05 
-1-0.20 



1.5 ±0.5 



17.95_o.o5 

J-0. ( l_o.49 
14 70+0-70 

< 13.87 

< 13.50 



Note: errors in the column densities correspond to a range of 
column densities; they are not the rms errors from fitting the 
Voigt profiles (see text). 
5(T upper limit. 



We measured ionic column densities in the three main 
components of the system and in a weaker one. The results 
are given in Table|K| The component where H2 is detected is 
not seen at all in Fell, but is prominent in the Cll* profile 
and detected in Sin, Sll and Znll (see Fig. |^. A careful 
modelling of the central blend leads to an accurate determi- 



nation of the column densities of the latter ions, and to an 
upper limit on the column density of Fe 11 in the component 
at Zabs ~ 2.59474 which is likely associated with the H2 
absorption lines (see Table |SJ . We compute overall metal- 
licities [Zn/H] = -1.02 ± 0.12 and [S/H] = -1.07 ± 0.10, 
and a mean depletion factor [Zn/Fe] = 0.31 ± 0.08. As pre- 
viously in the case of the DLA system toward Q 0347—383, 
the depletion is found to be larger in the component where 
H2 is detected ([Zn/Fe] > 0.70 at z^bs = 2.59474) relative 
to the other components ([Zn/Fe] = 0.38 ± 0.11, 0.25 ± 0.06 
and 0.28 ±0.05 at, respectively, Zabs = 2.59440, 2.59464 and 
2.59485). 



3.3 Q 1444±014, Zabs = 2.087 

Voigt-profile fitting to the Lyman-a line of this absorption 
system leads to a total neutral hydrogen column density 
log Ar(Hi) = 20.07 ±0.07 (see Fig.|7|l. Although this system 
would thus not qualify as a DLA system following the con- 
ventional definition, strong associated H2 lines are detected. 
Moreover, Aim, Siiv and Civ lines are barely detected at 
Zabs ~ 2.087 which shows that the amount of highly ionized 
gas is small in the system. The overall metal line profiles are 
unusual in displaying a series of well-separated sub-systems 
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Table 5. Ionic column densities in individual components of the 
DLA system at Zaba = 2.595 toward Q 0405-443 



Ion 


Transition 


log N ± o-iog jv 


b±at 




lines used 




(km s-i) 




2.59440 






Cii* 


1335.6,1335.7 


< 12.41 ± 0.14 


3.7 ±0.2 


Si 11 


1808 


14.35 ±0.05 


" 


Sii 


1250,1259 


13.91 ±0.04 




Fell 


2249,2260,2374 


13.98 ±0.02 




Zn II 


2026 


11.50 ±0.11 


" 


■^abs 


2.59464 






Cii* 


1335.6,1335.7 


13.09 ±0.10 


9.8 ±0.2 


Si II 


1808 


15 11 + 02 


Jl 


Sii 


1250,1259 


14.72 ± 0.02 


55 


Fell 


2249,2260,2374 


14.72 ±0.04 


55 


Znii 


2026 


12.11 ±0.04 


55 


■^abs — 


2.59474 






Ci 


1656 


< 11.84 " 




Ci* 


1656 


< 11.95 <= 




Ci** 


1657 


< 11.82 <= 




Cii* 


1335.6,1335.7 


13.20 ±0.20 


2.9; 


Si II 


1808 


14.16 ± 0.14 


55 


Sii 


1250,1259 


14.07 ± 0.09 


55 


Fell 


2249,2260,2374 


< 13.68 


55 


Znii 


2026 


11.52 ±0.13 




■^abs 


2.59485 






Ci 


1656 


^ 11.95 ± 0.34: 




Ci* 


1656 


< 11.95 




Ci** 


1657 


< 11.82 " 




Cii* 


1335.6,1335.7 


13.25 ±0.08 


7.2 ±0.3 


Si II 


1808 


15.18 ±0.02 




Sii 


1250,1259 


14.73 ±0.02 


55 


Fell 


2249,2260,2374 


14.75 ± 0.04 


55 


Zn II 


2026 


12.17 ±0.03 


55 



^ Also 1608, 1611. 

^ This component is likely associated with the H2 absorption 
lines observed at z^bs = 2.59471, 2.5 km s~^ blueward of 
2abs = 2.59474. 

5(7 upper limit. 

Blend. 

spread over ~ 400 km s~^. However, the Nl, Pll and Znll 
profiles suggest that most of the neutral hydrogen is con- 
centrated within ~ 40 km s~^ of the central clump (see 
Fig. |HJ . This is confirmed by the detection of several C I 
lines from different excitation states in two components at 
-Zabs = 2.08685 and 2.08697 (velocity separation 12 km s~^; 
see Fig.|UJ. The component at Zabs = 2.08697 is narrow and 
unresolved; the measured 6- value is 1.1 km s^^ (see TablelSJ. 

Molecular hydrogen lines from the J = 0, 1, 2 and 3 
rotational levels are detected in two different components 
most clearly seen in J ^ 2 (see Fig. HOI . Damping wings are 
definitively absent in J = and 1. The redshifts of the H2 
components are consistent with those of the C I lines. Since 
the latter are accurately determined, we simultaneously fit- 
ted all unblended H2 lines from different rotational levels 
with two components located at the measured redshifts of 
the C I lines. For a given component, the same broadening 
parameter was used independently of the J level of the lines. 
The column density in each J level was derived from several 



Table 6. Ionic column densities in individual components of the 
sub-DLA system at z^bs = 2.087 toward Q 1444+014 



Ion 


Transition 


log N ± CTlog iV 


b±at 




lines used 




(km s-i) 


■^abs 


2.08667 






Ni 


1134.1,1134.4,1134.9 


13.70 ± 0.10 


6.1 ± 0.9 


Alii 


1670 


11.77 ±0.07 




Si II 


1304,1526,1808 


13.55 ± 0.08 




Sii 


1253,1259 


13.79 ± 0.11 




Fell 


1121,1143,1608 


13.18 ± 0.09 




Znii 


2026,2062 


11.43 ± 0.18 




■^abs 


2 08679 






Ci 


1328,1560,1656 


12.67 ±0.09 


11.6 ± 2.8 


Ci* 


1656,1657.3,1657.9 


12.52 ± 0.14 


" 


Cii* 


1335.6,1335.7 


13.12 ±0.08 




Ni 


1134.1,1134.4,1134.9 


14.67 ±0.04 


6.9 ±0.6 


Mgii 


1239,1240 


< 14.61 


" 


Alii 


1670 


12.56 ± 0.03 




Si II 


1304,1526,1808 


14.15 ± 0.06 




Pii 


1152 


12.76 ± 0.06 




Sii 


1253,1259 


14.32 ± 0.05 




Till 


1910.6,1910.9 


< 12.10 ^ 




Crll 


2056,2062,2066 


12.14 ±0.12 




Fen 


1121,1143,1608 


13.83 ± 0.04 




Niii 


1317,1370 


12.31 ± 0.21 


„ 


Znll 


2026,2062 


11.66 ± 0.12 




■^abs — 


2.08692 






Ci 


1328,1560,1656 


12.82 ± 0.11 


1.1 ± 0.3 


Ci* 


1656,1657.3,1657.9 ^ 


12.42 ± 0.12 




Cii* 


1335.6,1335.7 


12.78 ± 0.20 




Ni 


1134.1,1134.4,1134.9 


14.46 ± 0.05 


5.1 ± 0.5 


Alii 


1670 


11.74 ±0.10 




Sin 


1304,1526,1808 


13.82 ±0.06 




Pii 


1152 


12.60 ±0.08 




Sii 


1253,1259 


14.17 ±0.06 




Fen 


1121,1143,1608 


13.23 ±0.09 




Zn n 


2026,2062 


11.77 ±0.07 





Note: the two components of the C I, C I * and C 11 * profiles are 
redshifted by ~ 5 km compared to what is observed for the 
other metal lines (whose redshifts are indicated). 

Also 1329.08, 1329.10, 1329.12, 1560.6, 1560.7. 

5cr upper limit. 
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Figure 7. Portion of the normalized UVES spectrum showing 
the damped Lya line of the sub-DLA system at z^bs = 2.087 
toward Q 1444+014. The best-fitted Voigt profile superimposed 
on the data corresponds to Af(Hl) = (1.2 ± 0.2) X lO^'^ cm"^. 
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Relative velocity (km s~^) Relative velocity (km s~^) 

Figure 8. Velocity profiles of lines from low ions and neutral species in the sub-DLA system at ^^bs = 2.087 toward Q 1444+014. 
Our best-fitting model of the central part of the profiles is superimposed on the spectra with vertical lines marking the location of 
individual components. The dotted parts in some of the synthetic profiles correspond to other transitions than the ones indicated. The 
two components where H2 is detected are labelled with number 1 and number 2 respectively. 



trials of Voigt-profile fitting to take into account the range 
of possible b values. Errors in tiie column densities therefore 
correspond to a range of column densities and not to the rms 
error from fitting the Voigt profiles. The results are given in 
Table 13 

The line profiles of S 11 and Zn 11 on the one hand, and 
Sill and Fell on the other, are different in the central clump 



(see Fig.|HJ. Whereas the S 11 and Znll profiles have a rectan- 
gular shape with two main components, the redder of these 
two components is barely detected in the wing of the clump 
in both SillA1304 and FenA1608. This demonstrates that 
Si and Fe depletions are much larger in the latter compo- 
nent. The redshifts of the main components of the system, in 
addition to a third one, were determined from simultaneous 
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Figure 9. Absorption lines from tlie ground-state of Cl (upper part of left panel) and fine-structure levels of Cl and C II (right panels) 
a-t 2abs = 2.087 toward Q 1444-1-014. The two fitted components, where H2 is detected, are indicated as in Fig. |H]by numbers in addition 
to vertical lines. The dotted parts in some of the synthetic profiles correspond to other transitions than the ones indicated. The origin 
of the velocity scale is the redshift of the main component as measured from the other metal lines (z^bs = 2.08679). The redward shift 
of the Cl components is apparent. The velocity profile of OlA1302 is shown for comparison. Note that AllllA1854 is barely detected 
(bottom part of left panel) in this sub-DLA system. 



fitting of most of tiie observed metal lines but independently 
of tiie fitting to the C I lines. This shows an apparent velocity 
offset of ~ 5 km s~^ between C I and other metal lines. The 
ionic column densities in the component in the red wing of 
the clump could also be derived accurately for a few species 
in the case of refractory elements using several transition 
lines covering a range of oscillator strengths (see Table [SJ. 
We find [Zn/Si] = 0.41 ± 0.13 and [Zn/Fe] = 0.69 ± 0.13 at 
Zabs = 2.08679, and [Zn/Si] = 0.85 ± 0.09 and [Zn/Fe] = 
1.40 ± 0.11 at Zabs = 2.08692. This is reminiscent of warm 
and cold Galactic disc cloud dust-depletion patterns respec- 
tively. The large depletion factors measured in the compo- 
nent at ZahB ~ 2.08692 are of the same order of magnitude 
as those measured in the component at 2abs = 1.96822 to- 
ward Q 0013-004 (Petitjean et al. 2002). This shows that, 
even though the depletion factors averaged over the entire 
profiles of DLA systems are usually moderate, they can be 
quite large in individual components. 



4 THE H2-SURVEY SAMPLE 

In this Section, we discuss all H2 measurements performed 
to date at high spectral resolution in DLA and sub-DLA 
systems at high redshift (zabs > 1-8). These measurements 
are summarized in Table |H| The H2 and Hi column densi- 



ties are the total ones integrated over the systems, except 
for Q 0013— 004 for which two sub-systems could be stud- 
ied individually. In most of the systems, there are several 
metal components and therefore the values for individual 
components may be scattered around the quoted mean. For 
each of the systems, we also give the Fe abundance, [Fe/H], 
and the metallicity, [X/H], relative to Solar abundances. The 
metallicity was measured from the abundance of Zn when 
Znll is detected, and we used either S or Si otherwise. As 
Zn is known to be little or undepleted onto dust grains, 
the [Zn/Fe] ratio should be a good indicator of dust de- 
pletion. S itself is undepleted onto dust grains. However, 
possible overabundance of a-elements relative to Fe at low 
metallicity may lead to overestimate the dust-depletion fac- 
tor when using [S/Fe]. However, the mean overabundance 
of a-elements should in general be less than about 0.25 dex 
(see Prochaska & Wolfe 2002). Si can be depleted onto dust 
grains and, hence, the ratio [Si/Fe] may lead to underesti- 
mate the dust-depletion factor, but probably by less than 0.3 
dex (see Petitjean et al. 2002). In the following, however, we 
do not attempt to derive exact values of the dust-depletion 
factor but only trends in the global behavior of the charac- 
teristics of the whole population. It will be shown that these 
considerations have little impact on the global results. 

In Figs. 111! and 1121 we compare the characteristics of 
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Figure 10. Voigt-profile fitting to the transition lines from the J = 0, 1, 2, 3 and 4 rotational levels of the vibrational ground-state 
Lyman band of H2 at z^bs ~ 2.087 toward Q 1444+014. H2 is detected on this line of sight in two gas clouds spanning 12 km in 
velocity space. 



Table 7. Voigt-profile fitting results for different rotational lev- 
els of the vibrational ground-state Lyman band of H2 toward 
Q 1444-f 014 



■^abs 


Level log A'' di cTiog AT 






(H2) 


(km s~^) 



2.08680 



2.08696 



15.68 
16.37 
15.31 
14.80 
<14 

< 13 
17.46 
18.03 
16.63 
15.26 
<:i4 

< 13 



-1-0.22 

0.08 
-1-0.27 
0.11 
-1-0.06 
0.03 
-t-0.06 
-0.04 
15 

75 ^ 

-1-0.10 
0.11 
-t-0.12 
0.12 
-t-0.57 
0.46 
-1-0.07 
0.08 
20 
75 ^ 



5.5 ±0.5 



2.6 ±0.6 



Note: errors in the column densities correspond to a range of 
column densities; they are not the rms errors from fitting the 
Voigt profiles (see text). 

Possible blends. 
^ 5cr upper limit. 



our sample of 33 systems (called in the following sample SI) 
with those of the global population of DLA systems. For 
the latter, we use the result of merging our sample with the 
sample of Prochaska et al. (2001). The corresponding sample 
will be called sample Spop and comprises 60 systems. In 
Fig. 111! different distributions for the global population of 
DLA/sub-DLA systems (sample Spop), sample SI and the 
sub-sample of systems where H2 is detected (sample SH2) 
are shown. 

The neutral hydrogen column density distributions are 
shown in the upper panel of Fig. 1111 There are 28 (resp. 32) 
systems with logA'^(Hl) smaller (resp. larger) than 20.6 in 
sample Spop, out of which 50 (resp. 59) percent belongs to 
sample SI. The two samples are statistically indistinguish- 
able although the mean H I column density is slightly smaller 
in our sample. The Kolmogorov-Smirnov test probability 
that the two distributions are drawn from the same parent 
population is Pks = 0.91 (two-sided case). It is also appar- 
ent, although the statistics are based on a smaller data set, 
that there is no systematic correlation between the detec- 
tion of H2 lines and the Hi column density. The log A'^(Hl) 
distributions are similar in sub-sample SH2 and in sample 
SI (Pks > 0.9999). 

It can be seen in the middle panel of Fig. 1111 that the 
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Table 8. Summary of molecular and metal contents in the H2-survey sample of DLA/sub-DLA systems 



QSO 






log7V(Hi) 


J 


log 
= 


Af(H2) 1 
J = 1 


log / 2 


[Fe/H] 3 


[X/H] 3 


X 


0000-263 


4.10 


3, 


.390 


21.41 + 0.08 ^ 


< 


13.9 


13.74:: ^ 


< -6.98 


-2.05 + 0.09 = 


-2.05 + 0.09 " 


Zn 


0010-002 


2.14 


2, 


.025 


20.80 ± 0.10 


< 


14.0 


< 14.2 


< -6.09 


-1.25 + 0.11 


-1.20 + 0.12 


Zn 


0013-004 


2.09 


1, 


.968 


^ 19.43 <^ 




16, 


77+0.05 d 
■ ' ' -0.07 


^ -2.36 


^ -2.33 


> -0.73 '^'^ 


Zn 


0013-004 


2.09 


1, 


.973 


20.83 + 0.05 d 




17, 


.72 / 20.00 <^ 


-2.81 / -0.64 d 


-1.75 + 0.05 <i 


-0.93 + 0.06 


Zn 


0058-292 


3.09 


2, 


.671 


21.00 + 0.10 


< 


13.8 


< 13.6 


< -6.69 


-1.76 + 0.10 


-1.42 + 0.11 


Zn 


0102-190 


3.04 


2, 


.370 


20.85 + 0.08 


< 


14.7 


< 14.6 


< -5.60 


-1.89 + 0.13 


-1.73 + 0.14 


Zn 


0112-306 


2.98 


2, 


.418 


20.37 + 0.08 


< 


13.8 


< 14.2 


< -5.72 


-2.50 + 0.09 


-2.31 + 0.08 


Si 


0112-306 


2.98 


2, 


.702 


20.15 + 0.07 


< 


14.1 


< 14.0 


< -5.50 


-0.89 + 0.10 


-0.33 + 0.11 


Si 


0112+029 


2.81 


2, 


.423 


20.70 + 0.10 


< 


13.6 


< 13.9 


< -6.32 


-1.35 ±0.11 


-1.14 + 0.13 


S 


0135-273 


3.21 


2, 


.800 


20.80 + 0.10 


< 


13.8 


< 13.8 


< -6.40 


-1.54 + 0.15 


-1.29 + 0.17 


S 


0347-383 


3.22 


3, 


.025 


20.56 + 0.05 




14, 


.55 + 0.09 


-5.71 + 0.10 


-1.72 + 0.06 


-0.98 + 0.09 


Zn 


0405-443 


3.02 


2, 


.550 


21.00 + 0.15 


< 


13.9 


< 13.7 


< -6.59 


-1.52 + 0.15 


-1.17 + 0.16 


Zn 


0405-443 
0405-443 


3.02 
3.02 


2, 
2, 


.595 
.621 


20.90 + 0.10 
20.25 + 0.10 


< 


18, 
13.5 


^g+0.21 
■^"-0.06 

< 13.5 


-2.44+0-?3 
—0.12 

< -6.15 


-1.33 + 0.11 
-2.15 + 0.10 


-1.02 + 0.12 
-1.83 + 0.10 


Zn 
Si 


0528-250 
0551-366 
0841 4-1 29 


2.78 
2.32 
2.50 


2, 
1, 
2, 


.811 
.962 
.374 


21.10 + 0.10 
20.50 + 0.08 ° 
2(1 QO + 1 


18, 
17, 
14.56:: 


99+0.23 5 
■^^-0.16 

.9+0.63 c 
■^^-0,90 
< 14.0 


-2.58+!|-?5 5 

—0.19 
9 70+0,64 c 
^' ' °— 0.90 

< —5.93 


-1.26 + 0.10 
-0.96 + 0.09 
— 1.71 + 0.11 


-0.75 + 0.10 
—0.13 + 0.09 


Zn 
Zn 
Zn 


0841 +1 29 


2.50 


2, 


.476 


20 65 + 1 


< 


14.0 


< 14.0 


< —6.05 


— 1.51 + 0.11 


— 1.52 + 0.11 


g 


1 037—270 


2.23 


2, 


.139 


iqvn + niof 


< 


13.7 f < 13.7 * 


< —5.40 ' 


„r) fjK + n 1 n f 


_n 26 + n 1 1 f 


Zn 


1101—264 


2.14 


1, 


.839 


1 9 35 + 04 


< 


14.0 


< 14.0 


< —4.75 


— 1 36 + 05 


—0 82 + 14 


g 


1117 1 

-L-L 1 / — loo 


o.yu 


3, 


.351 


90 ec; -1- n in s 


< 


13.7 


< 14.0 


— U.O 1 




— _L.ZO ^ U.IO 


Zn 


1157+014 


1.99 


1, 


.944 


21.70 + 0.10 


< 


14.3 


< 14.5 


< -6.69 


-1.73 + 0.10 


-1.32 + 0.10 


Zn 


1223+178 


2.94 


2, 


.466 


21.40 + 0.10 


< 


14.0 


< 14.0 


< -6.80 


-1.70 + 0.10 


-1.63 + 0.11 


Zn 


1232+082 


2.58 


2, 


.338 


20.90 + 0.10 




^ 17.19 ''•^ 


^ -3.41 


-1.73 + 0.13 


-1.21 + 0.15 


Si 


1337+113 


2.92 


2, 


.508 


19.95 + 0.05 


< 


13.8 


< 14.1 


< -5.37 


-2.05 + 0.07 


-1.59 + 0.09 


Si 


1337+113 


2.92 


2, 


.796 


20.85 + 0.08 


< 


13.6 


< 13.8 


< -6.54 


-2.02 + 0.09 


-1.69 + 0.11 


Si 


1444+014 


2.21 


2, 


.087 


20.07 + 0.07 




18.30^0-2^ 


_1 48+0-38 
J^-^O-0.38 


-1.58 + 0.09 


-0.60 + 0.15 


Zn 


1451+123 


3.25 


2, 


.469 


20.30 + 0.10 


< 


15.0: 


< 15.0: 


< -4.70 


-2.41 + 0.10 


-1.98 + 0.14 


Si 


1451+123 


3.25 


3, 


.171 


19.90 + 0.20 


< 


13.5 


< 13.5 


< -5.80 


-2.09 + 0.24 


-1.83 + 0.21 


Si 


2059-360 


3.09 


2, 


.508 


20.14 + 0.07 


< 


14.7 


< 14.5 


< -4.93 


-2.12 + 0.08 


-1.76 + 0.09 


S 


2059-360 


3.09 


3, 


.083 


20.85 + 0.08 


< 


13.5 


< 13.7 


< -6.64 


-1.83 + 0.09 


-1.65 + 0.12 


s 


2138-444 


3.17 


2, 


.852 


20.82 + 0.05 


< 


13.5 


< 13.1 


< -6.87 


-1.68 + 0.05 


-1.48 + 0.05 


Zn 


2332-094 


3.30 


3, 


.057 


20.30 + 0.08 


< 


12.8 


< 13.1 


< -6.72 


-1.47 + 0.08 


-1.49 + 0.20 


S 



^ We give total molecular hydrogen column densities summed up over all J levels in case of detection and upper limits for J = and 1 
in case of non-detection. 

2 Molecular fraction / = 2Ar{H2)/{2Af(H2) + N{Ul)). 

3 Fe abundances, [Fe/H], and metallicities, [X/H] with either X = Zn, or S, or Si. In order to derive meaningful abundance ratios, the 
components taken into account in the profiles are those simultaneously detected in the weak lines of Zn II, S II and Si II (see Ledoux et 
al. 2002a). 

■* [Zn/Fe] = 1.60 + 0.04. 

^ The fitting was performed on the new UVES data taking into account the two-component structure of the cloud. 

^ The total H2 column density of this system will be derived from higher quality data. The value derived by Srianand et al. (2000) 

should be considered as a lower limit. 

References: hydrogen column density and metal abundance measurements are from this work, unless otherwise indicated: (a) Lu et 
al. (1996); (b) Levshakov et al. (2001); (c) Molaro et al. (2001); (d) Petitjean et al. (2002); (e) Ledoux et al. (2002b); (f) Srianand & 
Petitjean (2001); (g) Peroux et al. (2002); (h) Srianand et al. (2000). 



metallicity distributions are similar in sample SI and in the 
global population Spop (Pks > 0.9999). There are 31 (resp. 
29) systems with [X/H] smaller (resp. larger) than —1.3 in 
sample Spop, out of which 55 (resp. 55) percent are part of 
sample SI. The mean metallicity is slightly smaller in our 
sample however. This could be a consequence of our sample 
including a large number of lines of sight with several DLA 
systems. Indeed, if these systems had a high metallicity and 
therefore a large dust content, they would have obscured 
the background quasars. However, this does not mean that 
these lines of sight are biased against the presence of H2 . We 
indeed detect H2 at Zabs ~ 2.595 toward Q 0405—443, which 



is a line of sight containing no less than three DLA systems. 
The important result is that the [X/H] distributions are dif- 
ferent in sub-sample SH2 and in sample SI (Pks = 0.03): 
systems where H2 is detected are definitively amongst the 
most metal-rich. Interestingly, none of the 17 systems with 
[X/H] < —1.3 in sample SI show detectable amount of H2. 
However, H2 is detected down to [X/H] ~ —1.2. 

The distributions of depletion factors (Fig. 1111 lower 
panel) are also not significantly different in sample SI com- 
pared to sample Spop (Pks = 0.98). There are 35 (resp. 
25) systems with [X/Fe] smaller (resp. larger) than 0.4 in 
sample Spop, out of which 60 (resp. 48) percent belongs to 
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Figure 11. Distributions of neutral hydrogen column densities 
{upper panel), metallicities {middle panel, measured from either 
[Zn/H], or [S/H], or [Si/H]) and depletion factors (Jower panel; 
measured from either [Zn/Fe], or [S/Fe], or [Si/Fe]) in, respec- 
tively, the global population of DLA/sub-DLA systems (overall 
distribution), our H2-survey sample (dotted histogram) and the 
sub-sample of H2-detected systems (hashed histogram). 
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Figure 12. Neutral hydrogen column density {upper plot) and 
depletion factor (iowerpJot) versus metallicity in, respectively, the 
global population of DLA/sub-DLA systems (all symbols), our 
H2-survey sample (all squares) and the sub-sample of H2-detected 
systems (dark squares). The typical error bars are shown in the 
upper-left part of each plot. A correlation between metallicity and 
depletion factor is present at the 4a significance level, confirming 
the trend previously detected by Ledoux et al. (2002a). 



sample SI. Although H2 is detected in a DLA system with 
[Zn/Fe] as low as 0.3 (at z^bs = 2.595 toward Q 0405-443), 
there is a clear tendency for depletion factors to be larger in 
systems where H2 is detected. The [X/Fe] distributions are 
different in sub-sample SH2 and in sample SI (Pks ~ 0.02). 
In particular, H2 is detected in the five systems with the 
largest depletion factors ([X/Fe] ^ 0.7). 

In Fig. 1121 the H I column density (upper panel) and the 
depletion factor (lower panel) are both plotted versus metal- 
licity. Measurements from sample SI are shown by squares, 
with dark squares for systems where H2 is detected. It is 
apparent that, in the overall population of DLA systems, 
there is a lack of systems with both a high metallicity and 
a large H I column density. This was first noticed by Boisse 
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Figure 13. Mean H2 molecular fraction / = 2Af(H2)/(27V(H2) + A'^{Hl)) versus neutral hydrogen column density. Measurements in 
DLA/sub-DLA systems are indicated by dark squares for H2 detections and shaded ones for upper limits. Observations along lines of 
sight in the Galaxy (Savage et al. 1977) and the LMC and SMC (Tumlinson et al. 2002) are indicated by, respectively, asterisks, circles 
and triangles. 



et al. (1998) who suggested that this might be due to the 
fact that, if these systems exist, they may avoid detection 
because the background quasars are attenuated by dust ex- 
tinction (Fall & Pei 1993). This assumes at least some corre- 
lation between metallicity and the amount of dust. It can be 
seen in the lower panel of Fig. ll2l that indeed there is a trend 
for the most metal-rich systems to have larger depletion fac- 
tors. This is detected at the 4a significance level using a 
Kendall rank correlation test on censored data and taking 
into account the measurement uncertainties. This trend was 
already noticed before by Ledoux et al. (2002a), at a lower 
significance level, probably due to a less homogenous sample, 
larger measurement uncertainties and the observation of a 
narrower metallicity range. However, from this alone, it was 
unclear whether this trend is a consequence of dust depletion 
effects or peculiar nucleosynthesis history of zinc (Prochaska 
& Wolfe 2002). However, a correlation is observed in both 
our sample and also the global population of DLA systems. 
Moreover, it is clear from Fig. 1121 that H2 is detected at 
both the highest metallicities and the largest depletion fac- 
tors. This strongly favors the fact that the correlation be- 
tween metallicity and depletion factor is a consequence of 
dust depletion effects. 



5 DISCUSSION 

5.1 Molecular hydrogen in DLA systems 

If we consider only DLA systems, with logA'^(Hl) ^ 20.3, 
firm detection of H2 is achieved in 6 out of 24 systems. If we 
also include sub-DLA systems, this amounts to 8 detections 
out of 33 systems. Therefore, H2 is detected in more than 20 
percent of the DLA systems. However, three of the systems 



where H2 is seen were already known. If we exclude these sys- 
tems, we actually detect H2 in 13 to 20 percent of the newly 
surveyed systems. Note that we do not include in these cal- 
culations the tentative detections toward Q 0000—263 (Lev- 
shakov et al. 2000) and Q 08414-129 (Petitjean et al. 2000). 

The detection probability of H2 with column densities 
greater than 10^* cm~^ is more than 90 percent along lines of 
sight through the ISM of our Galaxy and in the SMC (Sav- 
age et al. 1977; Richter et al. 2001; Tumlinson et al. 2002). 
However, only 52 percent of the lines of sight through the 
LMC have detectable H2 fines (Tumlinson et al. 2002). These 
differences in the detection probability cannot be ascribed 
solely to different chemical enrichment histories as the ISM 
of the LMC has a mean metallicity only 0.3 dex lower than 
the Galactic local ISM and the metallicity of the SMC is 
0.6 dex lower than that (Russell & Dopita 1992). Therefore, 
the small detection probability in DLA systems is probably 
not only related to their low metallicities. The level of star- 
formation activity in or nearby the systems might well be 
the primary reason for the low detection rate of H2 as it has 
been advocated in the case of the LMC (Kim et al. 1999). 

As noted above, in our sample the detection of H2 is 
not correlated with the H I column density. The H I column 
density distribution is similar in sample SI and in the sub- 
sample of systems where H2 is detected (see Fig. IllH . In 
particular, H2 is detected in two sub-DLA systems, with 
logA'^(Hl) < 20.3. In addition, amongst the systems where 
H2 is detected there is no correlation between the logarithm 
of the molecular fraction, log/, and logA'^(Hl) (see Fig. ll3ll . 

In Fig. 1131 confirmed detections and upper limits on the 
molecular fraction in sample SI are represented by, respec- 
tively, dark and shaded squares. Observations along lines of 
sight in the LMC and SMC from Tumlinson et al. (2002) 
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Figure 14. Mean H2 molecular fraction versus depletion onto 
dust grains as estimated from the abundance of either X=Zn, or 
S, or Si relative to Fe, in the H2-survey sample (see Table IS! . H2 
detections, measured by our group, are indicated by dark squares. 
The dotted line segment links data points for the DLA system 
toward Q 0551— 366 considering either [Fe/Zn] or [Fe/S]. This is 
indeed the only one case where the difference between measured 
[X/Fe] ratios is significantly larger than 0.1 dex (i.e. 0.24 dex). 
The two dotted lines are drawn for illustrative purposes (see text). 




Figure 15. H2 molecular fraction versus the amount of dust, or 
dust-to-gas ratio kx = 101^/^1(1 _ ;^q[Fc/x]-)^ ^^^^i either X=Zn, 
or S, or Si. Drawing conventions are the same as in Fig. 1141 A weak 
trend between the two quantities is apparent; however, several 
orders of magnitude spread in /(H2) for a given dust content, 
among the detected cases, suggests that in addition to the amount 
of dust the physical conditions of the gas (density, temperature, 
UV flux) play an important role in governing the formation of H2 
in DLA systems. 



and in the Galaxy from Savage et al. (1977) are also shown. 
It is apparent that most DLA systems are similar to the 
LMC and/or SMC lines of sight. In particular, and con- 
trary to what happens in our Galaxy, a significant num- 
ber of lines of sight through the LMC have log / < —6 
for 20.5 ^ logAr(Hl) ^ 21.5. This is also the case for 58 
percent of the DLA systems (20.3 < logiV(Hl) 21.7) in 
which we could achieve such a limit, which means that this 
number is a lower limit on the fraction of systems having 
such a small molecular fraction. Indeed, if all detection lim- 
its are considered to be that constraining then the fraction 
of systems with log / < —6 is 75 percent. Note that several 
systems have a very small molecular fraction, log / < — 7. 
Interestingly, this is an order of magnitude smaller than the 
primordial freeze-out fraction of H2 molecules (Lepp & ShuU 
1984). It is very difficult to reproduce such small values 
in models of cool clouds (see Petitjean et al. 2000; Liszt 
2002). When H2 is detected, the molecular fraction is in the 
range —3 < log/ < —1 with some scatter, with the ex- 
ception of the DLA system toward Q 0347—383 for which 
log / = —5.71. The former values are similar to what is seen 
along the LMC and SMC lines of sight. 

The sharp transition from log / < — 4 to log / > — 2 at 
a column density logA*'(IIl) ~ 20.7 observed for the Galac- 
tic lines of sight (Savage et al. 1977) is not clearly present 
in the SMC, in the LMC nor in sample SH2. This transi- 
tion is believed to occur at the point at which the clouds 
achieve self-shielding. However, a characteristic H I column 
density can possibly correspond to some self-shielding scale 
only if the parameters controlling the formation equilibrium 



of H2 are similar along different lines of sight, with very little 
scatter, which may be the case in our Galaxy. Therefore, our 
results suggest that DLA systems span a wide range of phys- 
ical conditions. This is also suggested by the large scatter 
seen in the metallicity measured in DLA systems at a given 
redshift (see e.g. Prochaska & Wolfe 2002). This conclusion 
is strengthened by the fact that there is a clear dichotomy 
in log / values whatever the H I column density is, with very 
few measurements in the range — 5 < log / < — 3. This bi- 
modal distribution could be a consequence of self-shielding, 
which would mean that, in DLA systems, self-shielding can 
be achieved at any H I column density probably because the 
UV radiation field varies strongly from one system to the 
other. Alternatively, as shown by Tumlinson et al. (2002) 
in the case of the Magellanic Clouds, most of the observa- 
tions can be explained by models where the formation rate 
of II2 onto dust grains is reduced and the ionizing flux is 
enhanced relative to what is observed in our Galaxy (see 
also the models by Shaye 2001). Basically, the above char- 
acteristic Hi column density, logA*'(Hl) — 20.7, observed 
for Galactic clouds, would be located beyond the observed 
range in the case of DLA systems (i.e. log A'^(Hi) > 21.8). 

Note that there are two sub-DLA systems, with 
log7V(Hi) < 20.3, having detected H2 molecules. These sys- 
tems, at Zabs ~ 1.968 toward Q 0013— 004 and Zabs = 2.087 
toward Q 1444-1-014, are amongst the most metal (and dust) 
-rich absorbers of sample SI. They, together with the DLA 
system at Zabs = 2.338 toward Q 1232-1-082, have the high- 
est molecular fractions, and could be part of the population 
arising in disc-like gas. 



18 Cedric Ledoux, Patrick Petitjean, & R. Srianand 



5.2 The role of dust 

It can be seen on Fig. 1111 that the probabihty of detect- 
ing H2 is unity for systems with large depletion factors, 
[X/Fe] ^ 0.7. This is a crude indication that the presence 
of dust is an important factor for the formation of H2 in 
DLA systems. It is known that H2 formation is mediated by 
either dust grains if the gas is cool and dense, or H~ if it 
is warm and dust-free. If the former process is dominant in 
DLA systems then some correlation between H2 molecular 
fraction and amount of dust is expected. However, in the 
case of DLA systems it is difficult to quantify the amount 
of dust that is present. We first plot in Fig. 1141 molecular 
fractions versus depletions (probably onto dust grains) as 
estimated from the abundances of Zn relative to Fe, or S 
when Zn is not detected, or Si when Zn is not detected and 
sulfur lines are blended. A trend between both quantities is 
present at the 3.6(j significance level using a Kendall rank 
correlation test taking into account upper and lower limits 
on log /. Note however that there is a two orders of magni- 
tude scatter in log / for a given depletion factor so that there 
is basically no correlation between log / and [X/Fe] if only 
detections are considered (correlation significant at the < lo 
level only). The above trend is therefore a consequence of 
H2 being detected in none of the systems with [X/Fe] < 0.4 
except one (at z^ha = 2.595 toward Q 0405— 443). The corre- 
sponding systems have a very low dust content. A depletion 
factor [Zn/Fe] ~ 0.4 corresponds to a very small E{B — V) 
(~ 0.01) for physical conditions similar to the ones prevail- 
ing in our Galaxy. We have seen that the depletion factor 
is correlated with the metallicity in DLA systems. There is 
therefore some correlation between the depletion factor and 
the amount of dust. However, the scatter is large. This is why 
it can be, in general, misleading to use only the depletion 
factor as an indicator of the dust content of the absorbing 
gas. 

A better indicator of the dust content is the dust-to- 
gas ratio: Kx = 10''^/"'(1 - lO''''^''^'), where X stands for a 
reference element which is little or unaffected by dust de- 
pletion effects. In Fig. 1151 we plot the molecular fraction 
versus the amount of dust as measured by the dust-to-gas 
ratio. It can be seen that, as expected, low dust content, 
logAtx < -1.5 (e.g. [X/H] < -1.3 and [X/Fe] < 0.5), implies 
low molecular fraction, log / < —5, while larger molecular 
fraction, log / > —4, is only found for higher dust content, 
logKx > —1.5. However, large amounts of dust do not al- 
ways imply large molecular fractions. There are a few non- 
detections with log Kx > —1.5 and log / < — 5. In particular, 
the two systems with log kx — —0.5 and log / < —5 are sub- 
DLA systems (see also Fig. 1121 upper plot, at [X/H] ~ —0.3). 
Therefore, in these systems the lack of molecules could be 
the consequence of enhanced radiation field below the Ly- 
man limit. 

Another important factor governing the presence of 
H2 molecules appears to be the local physical conditions 
of the absorbing gas. Relatively large molecular fractions, 
log / > —4, are found in gas having large particle densities, 
riH > 20 cm~'^, and low temperatures, T < 300 K (see Pe- 
titjean et al. 2002; Ledoux et al. 2002b). One should keep 
in mind that the formation rate of H2 goes linearly with the 
density of dust grains while it goes as the second power of 
the H I density. Therefore, it is natural that, even though the 
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Figure 16. Upper panel: redshift distribution of the DLA/sub- 
DLA systems in our sample. The sub-sample of H2-detected sys- 
tems is indicated by the hashed histogram. Lower panel: H2 
molecular fraction versus redshift. Detections are shown by dark 
squares and upper Hmits by shaded ones. Although the statistics 
are based on small numbers, the probability of finding molecules 
S't 2^abs ~ 2 is highest. This probably reveals that the amount of 
dust in DLA systems increases with time. 



presence of dust is an important factor for the formation of 
molecules, local physical conditions such as gas density, tem- 
perature and local UV field play a major role in governing 
the molecular fraction of a given cloud in DLA systems. 

Note that the above discussion is based on the determi- 
nation of mean metallicities in DLA systems. In most cases, 
this means that the metallicity is calculated over some veloc- 
ity range and, therefore, over several components. However, 
we have shown that the depletion in the components where 
H2 is detected can be much larger than the mean depletion. 
As an example, a large depletion factor has been observed 
in a well-defined, weak H2 component in the DLA system 
at Zabs = 1.97 toward Q 0013-004 (Petitjean et al. 2002). 
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This kind of component can usually be revealed by the cor- 
responding, narrow C I absorption line. However, for other 
species the component is most often blended with other com- 
ponents in the rest of the profile. This means that absorption 
lines due to cold gas, i.e. related to C l components, can be 
hidden in the absorption line profile produced by a more per- 
vasive medium. This implies that absolute metallicities are 
not accurately determined in the cold gas. This fact is proba- 
bly not important for volatile elements as for them the dom- 
inant components actually produce the strongest absorption 
lines, but could be of importance for Fe co-production ele- 
ments which can be depleted onto dust grains. Consequently, 
absorption lines from the latter elements can be weak and 
lost in the profiles. Therefore, we could underestimate the 
depletion factors in these components. 

5.3 Evolution with redshift 

In the previous Section, we have given two new arguments in 
favor of the presence of dust in DLA systems. First, there is 
a correlation between metallicity and depletion factor and, 
second, molecular hydrogen is observed in the systems hav- 
ing the highest metallicities and the largest depletion fac- 
tors. In a simple picture of star-formation history, one can 
imagine that the gas is slowly enriched by on-going star for- 
mation and, therefore, metallicity and dust content should 
increase together with time. In Fig. 1161 we plot the to- 
tal number of DLA/sub-DLA systems in our sample versus 
redshift along with the number of systems where H2 is de- 
tected. Although the statistics are based on small numbers, 
the probability of finding molecules in the Zabs ~ 2 systems 
is highest. This probably indicates that the amount of dust 
in DLA systems increases with time. Although the increase 
of Zn metallicity with decreasing redshift is not fully estab- 
lished (e.g. Pettini et al. 1999), there is some evidence for 
such an increase at ^abs < 3.5 (Prochaska & Wolfe 2002; 
Kulkarni & Fall 2002). Consequently, the possible lack of 
increase with time of the mean Fe metallicity in DLA sys- 
tems (e.g. Prochaska & Wolfe 2002) could be related to the 
fact that the intrinsic increase in metallicity is hidden by 
larger depletion onto dust grains at low and intermediate 
redshifts (see Ledoux et al. 2002a). 



6 CONCLUSIONS 

Although the presence of dust in DLA systems was claimed 
very early (Pei, Fall & Bechtold 1991), the issue of whether 
the abundance pattern observed in DLA systems refiects 
depletion of refractory elements onto dust grains or nucle- 
osynthesis effects has remained controversial (see e.g. Lu 
et al. 1996; Pettini et al. 1997). Recently, several studies 
have shown that the two effects, dust depletion and pecu- 
liar nucleosynthesis history, should be invoked altogether to 
explain the observations (Vladilo 1998; Prochaska & Wolfe 
2002; Ledoux et al. 2002a). In any case, all these studies con- 
clude that the dust content of DLA systems is usually small. 
However, it is possible that the current sample of DLA sys- 
tems is biased against high-metallicity and dusty absorbers. 
Indeed, Boisse et al. (1998) have noticed that there is a lack 
of absorbers with both a large N(R l) and a high metallicity. 
If this dust-induced bias exists, however, it probably cannot 



lead to underestimate the H l mass in DLA systems by more 
than a factor of two (Ellison et al. 2001). 

In the course of a survey of 24 DLA and 9 sub-DLA sys- 
tems, we have confirmed four previous detections of molecu- 
lar hydrogen and made three new ones. The bulk of the DLA 
population is found to have low metallicities, [X/H] < —1.0, 
together with small depletion factors, [X/Fe] < 0.5. Al- 
though molecules can be found along lines of sight with small 
depletion factors ([Zn/Fe] = 0.3), they are generally found 
in the most metal-rich systems, [X/H] > —1.0, with large de- 
pletion factors, [X/Fe] > 0.5. This clearly demonstrates the 
existence of dust in at least 20 percent of the DLA systems. 
Moreover, we have found very large depletion factors in two 
sub-systems where H2 is detected, at Zabs = 1.96822 toward 
Q 0013-004 (Petitjean et al. 2002), where [Zn/Fe] = 1.6, 
and at Zabs = 2.08692 toward Q 1444-1-014 (see Sect.^SJ, 
where [Zn/Fe] = 1.4. This shows that dust depletion can 
be quite large in some of the DLA systems. The extinction 
is small in these cases however, probably because of rela- 
tively small Hi column densities (e.g. logA''(Hl) < 19.5 in 
the Zabs = 1.96822 system toward Q 0013-004). In addition, 
there is a correlation between metallicity and depletion fac- 
tor in all DLA samples (large depletion factors are seen in 
systems with high metallicities). This trend has been noticed 
before by Ledoux et al. (2002a). However, it was unclear 
whether this trend is really due to dust depletion. Prochaska 
& Wolfe (2002) indeed argued that such a correlation may 
arise from a significant dispersion in the production of zinc. 
However, the correlation is present in both our sample and 
also the global population of DLA systems. Moreover, it 
is clear that molecular hydrogen is mostly detected at the 
highest metallicities and the largest depletion factors. This 
strongly favors the hypothesis that the correlation between 
metallicity and depletion factor is actually a consequence 
of dust depletion effects. Not only dust could be present in 
most of the DLA systems, but also the possibility that the 
population of DLA systems is biased against the presence 
of systems having both a high metallicity and a large dust 
content still needs to be considered seriously. 

One of the striking results of the survey is that for most 
of the DLA systems the molecular fraction is very small 
(log / < —6). This is much smaller than what is observed in 
the disc of the Galaxy at similar Hi column densities (see 
Savage et al. 1977) . However, the situation is comparable to 
what is observed along lines of sight in the LMC and SMC 
(Tumlinson et al. 2002). Indeed, similar upper limits on the 
molecular fraction have been measured for a large number 
of lines of sight in the LMC. Moreover, the mean molecular 
fraction in DLA systems where H2 is detected is of the order 
of 0.01, about the same as in the LMC and SMC and about 
a factor of ten smaller than in the disc of the Galaxy. This 
is probably a consequence of reduced formation rate of H2 
onto dust grains because the temperature of the gas is high 
(Petitjean et al. 2000). Indeed, in the framework of models 
of cool clouds it is difficult to explain the very low values 
of the upper limits (log / < —7) achieved for a handful of 
DLA systems in our sample (see fig. 7 of Liszt 2002). This 
can also be due to an intense ambient UV fiux (see models 
by Tumlinson et al. 2002), but this is unclear and should be 
investigated in detail by analysis of the systems where H2 
is detected (see Srianand & Petitjean 1998; Petitjean et al. 
2002; Ledoux et al. 2002b). 
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Star-formation activity is probably not intense in the 
close vicinity of most of the systems. M0ller et al. (2002) 
found that a few DLA systems are associated with Lyman- 
break galaxies. However, from spectroscopy of 11 out of 22 
candidates in 6 fields, targeting 8 DLA systems, they iden- 
tified only one galaxy counterpart. This means that, even if 
DLA systems can be associated with Lyman-break galaxies, 
they must have most of the time moderate star-formation 
activity with luminosities L < L* (see e.g. Fynbo, M0ller 
& Thomson 2001). In addition, in the framework of current 
models of structure formation, DLA systems preferentially 
sample the outer regions of galaxies at the faint end of the 
galaxy luminosity function (Haehnolt ot al. 2000). One way 
of explaining the small molecular fractions observed in DLA 
systems is that the associated emitting object is faint and 
has moderate star-formation activity, but that the gas giv- 
ing rise to the DLA absorption line is located inside the 
regions where star formation occurs. The probably large tur- 
bulence induced by on-going star formation may then help 
explain the required large covering factors associated with 
such regions. This would imply that star formation is gentle 
and diffuse in the systems. The diffuse UV flux can be esti- 
mated from the excitation of molecular hydrogen. This will 
be done in a forthcoming paper. Additional deep imaging in 
the vicinity of the QSO lines of sight is also needed to unveil 
the nature of star-formation activity in DLA systems. 
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